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ABSTRACT 

In the dynamic diffusion limit of radiation hydrodynamics, advection dominates diffusion; the latter 
primarily affects small scales and has negligible impact on the large scale flow. The radiation can 
thus be accurately regarded as an ideal fluid, i.e., radiative diffusion can be neglected along with 
other forms of dissipation. This viewpoint is applied here to an analysis of simple waves in an ideal 
radiating fluid. It is shown that much of the hydrodynamic analysis carries over by simply replacing 
the material sound speed, pressure and index with the values appropriate for a radiating fluid. A 
complete analysis is performed for a centered rarefaction wave, and expressions are provided for 
the Riemann invariants and characteristic curves of the one-dimensional system of equations. The 
analytical solution is checked for consistency against a finite difference numerical integration, and the 
validity of neglecting the diffusion operator is demonstrated. An interesting physical result is that 
for a material component with a large number of internal degrees of freedom and an internal energy 
greater than that of the radiation, the sound speed increases as the fluid is rarefied. These solutions 
are an excellent test for radiation hydrodynamic codes operating in the dynamic diffusion regime. 
The general approach may be useful in the development of Godunov numerical schemes for radiation 
hydrodynamics . 

Subject headings: hydrodynamics, radiative transfer 



1. INTRODUCTION 

The equations of radiation hydrodynamics (RHD) are 
used to model fluids in which matter and radiation are 
coupled. In an opaque material, the coupling is strong 
and acts primarily to maintain the material and radia- 
tion fluids in local thermodynamic equilibrium (LTE). 
Once at a common temperature, the radiation acts as 
an additional source of pressure as well as a diffusive en- 
ergy sink, and the material advects the radiation. The 
relative importance of diffusion and advection depends 
upon both the optical depth of the system and its char- 
acteristic velocity. If (3t -C 1, where v = (3c is the fluid 
velocity and / = tA p is a characteristic length (A p being a 
photon mean- free path), diffusion dominates advection; 
such a system is referred to as being in the static diffu- 
sion limit. If (3t 1, advection dominates diffusion and 
the system is in the dynamic diffusion limit. This paper 
is concerned primarily with dynamic diffusion. 

In the dynamic diffusion limit, diffusion only operates 
on length scales A <C I, in a manner analogous to vis- 
cosity and heat conduction (albeit on different scales). 
The bulk flow properties can thus be modeled accurately 
by treating the radiation fluid as ideal, i.e., by neglect- 
ing radiative diffusion along with other forms of dissipa- 
tion. Many of the methods developed for analyzing fluids 
without radiation can thus be applied to a radiating fluid 
operating in the dynamic diffusion limit. The effects of 
diffusion are limited to boundary layers and shocks and 
other regions of the flow in which gradients are large. 

This viewpoint is pursued here by analyzing the equa- 
tions of ideal RHD in one dimension. I begin in ^ by 
outlining the equations of ideal RHD in the gray diffu- 
sion limit in a frame comoving with the fluid. A review 
of the thermodynamics of a radiating fluid is provided 
and the impact of radiation on a centered rarefac- 



ence numerical integration are included as a consistency 
check, along with an example of the effects of dissipa- 
tion. A discussion of how to generalize this approach to 
include more general boundary conditions is given in fj5l 
along with expressions for the characteristic curves and 
Riemann invariants of the combined fluid. A summary 
and a discussion of how these ideas might be used in the 
development of Godunov schemes for numerical RHD is 
given in Sj6] 

2. EQUATIONS 

The energy equations for an optically-thick gray 
medium in LTE, in a frame comoving with the fluid and 
neglecting scattering, are 
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where p, U and T are the mass density, internal energy 
density and temperature of the material, E is the energy 
density of the radiation, c is the speed of light, a is the 
radiation constant, k ~ A p 1 is the absorption opacity in 
units of inverse length, and d/dt is the derivative follow- 
ing a fluid clement. Equation |T]) is specific to an ideal 
gas equation of state with material pressure 



p m = (7 - m 



(3) 



where 7 is the adiabatic index of the material. 

The three relevant time scales for these equations are 
the advection time scale t a dv — l/v, the diffusion time 
scale tdiff = k1 2 /c and the time scale for coupling be- 
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mean time between photon collisions) 
scalings are 
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so that in the dynamic diffusion limit (r ^> f3 1 ) 

< ■C t dl ff. (7) 

Thermal equilibrium is thus established on a short time 
scale, and diffusion can be neglected for spatial variations 
A ~ I. This can be seen explicitly by subtracting equa- 
tion J2J from equation (JTJ) : 
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The diffusion operator is 0([/3t] _1 ) relative to the ad- 
vection terms, which are in turn 0(/3r _1 ) relative to the 
left hand side, so that 

E ~ aT 4 (9) 

in the dynamic diffusion limit. Expression will gener- 
ally be considered to be an equality throughout this anal- 
ysis, although of course departures from thermal balance, 
no matter how small, are necessary to maintain coupling 
between the radiation and the material. Denoting the 
total pressure P m + P r by P, where 
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E ^ aT 4 
~3 ~ IT' 
the equations of ideal RHD are given by 
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Coupled with the equations of state Q and (JXDJ) for the 
material and radiation, these equations form a closed set. 
Their form in one Cartesian dimension is the basis for the 
analysis in the following sections. 

3. THERMODYNAMICS 

Most of the contents of this section are not new; the re- 
sults can be found in standard textbooks (|Chandrasekharl 
[19671 ICox fc Giuhl [l968t [Mihalas fc Mihalad Il984[) and 
are only included here for completeness and ease of ref- 
erence, and to highlight the most general features of the 
solutions to be described in the next section. The sum 



of equations p3|) and (jl4|) is the first law of thermody- 
namics for the combined material-radiation fluid: 

d(U + E)- Uu+^E^ ^ =0. (15) 

It is straightforward to demonstrate that this is equiva- 
lent to 

d(s m + Sr )=0 1 (16) 

where 



(17) 



is (to within a constant) the specific entropy of the ma- 
terial and 

4aT 3 , . 

is the specific entropy of the radiation (i.e, the entropy of 
radiation per unit mass of material). 1 Changes in ideal 
RHD are adiabatic: the total entropy of the combined 
fluid is conserved for a fluid element. 

The first law combined with the equation of state for 
the material (U ~ pT) implies 
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(20) 



is the ratio of radiation to material pressure. The density 
can be expressed as a function of a by combining (Q3 
and ([201): 
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where a subscript denotes a reference value; for the prob- 
lem analyzed in the following section, it is the quiescent 
value before the rarefaction wave is excited. An impor- 
tant implication of the above expression is that rarefac- 
tion (compression) is associated with an increase (de- 
crease) in the pressure ratio a for 7 < 4/3 (i.e., p ^ po 
when a ^ ao); for 7 > 4/3, the opposite conditions hold. 
The other fluid quantities are given by 

m \ 37 — 4 / \ 7— 1 

— ] =(—) e 4(T- 1 )("-«o>, (23) 
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1 Notice that s r is the entropy of a photon gas coupled to a 
material fluid; it is not in general equivalent to s m with 7 = 4/3. 
The two are equivalent (to within a constant factor) only in the 
high energy density limit. 
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Expressions Jp2l) - ([25f match expressions 9.127 of 
ICox fc Ciulil (1968) for 7 = 5/3 and a = Z. 

For 7 = 4/3, a = const, and expressions ([2"2"|) - ([2"5|l 
must be replaced with 
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Small changes in the total pressure satisfy the following 
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expression: 
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Combining the above expression with expression |T 
gives the adiabatic sound speed for the radiating fluid 
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where the total entropy is held constant and 
d\nP 1 /_ [ 7 -l][l + 4a] 2 
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d\np l + a\* l + 12[7-l]a 
The closed form expression for C a cl function of a is 



c a = c I — 



2(3-y-4) 
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where To is Ti evaluated at a = uq. The above expres- 
sion can be combined with expression (|A2[) to obtain a 
closed form expression for P{p); one can thus regard the 
radiation in ideal RHD as a modification of the equation 
of state for the fluid. 

4. CENTERED RAREFACTION WAVE 

This section will attempt to follow closely the standard 
hydrodynamic analysis, a cogent expression of whi ch can 
be found in §§99-105 of lLandau fc Lifshitd (fl987h . One 
of the simplest one-dimensional flows to analyze is a cen- 
tered rarefaction wave. The prototypical example of such 
a flow is that generated behind a piston moving out of 
a semi-infinite fluid at constant velocity. Since there are 
no characteristic time or length scales in such a system, 
all of the flow variables must depend upon the coordi- 
nates through the similarity variable £ = x/t, so that 
d/dt = -{{/t)d/d£ and d/dx = {l/t)d/dt The form of 
equations (|11[) and p2[) for a one-dimensional similarity 
flow is 

. s dp dv , . 

(«- 0^ + ^=0, (33) 



dv 1 



= 0. 



(34) 



Combining these equations with equation (|30|) gives 
W ± Ca = £ (35) 

and 

v = ± ,h*£ = ±[*L, (36) 
P 7 pc a 



where a rarefaction wave corresponds to the plus sign and 
a compression wave to the minus sign. These correspond 
to the hydrodynamic results with c m = \J^P m j p — > c a 
and P m — > P. Details on evaluating the above integral 
can be found in the Appendix; an approximate solution 
valid for 7 > 1.1 can be obtained from the following 
considerations. 
From (|3TJj) and p)]) . 
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The above expression is exact. Since Ti varies much more 
slowly with ^ than either v or c a (7 < Ti < 4/3 and from 
|35| one expects v and c a to vary approximately linearly 
with £), it can be approximated by 
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where the velocity is defined to be zero in the initial 
state. This corresponds to the hydrodynamic result with 
c m ~ * c a and 7 — > T a . The velocity slope varies between 



2 dv 6 
7 + 1 < ^ < 7" 



(41) 



To obtain the solution for the other flow variables, it is 
necessary to solve expressions (|3"2"|) and P0)l implicitly to 
obtain an d then insert the result into expressions 

The wave front is located at x = coi and propagates 
away from the piston. A region of constant velocity equal 
to the velocity of the piston is located between the piston 
and the point 



Xfr 



Co 



r + 1 , 



(42) 



where v p < is the piston velocity. The location of 
this transition between the similarity flow and the re- 
gion adjacent to the piston travels away from the piston 
for \v p \ < 2co/(r + 1); otherwise it travels towards the 
piston. For a piston velocity greater in magnitude than 



p I crit 



r„ 



1 



c , 



(43) 



the fluid is evacuated between the piston and the point 
x = — \vp\ cr i t t (the fluid quantities go to zero). 

Figure [T] shows the slope of the velocity as a func- 
tion of ao, the quiescent ratio of radiation to mate- 
rial pressure, for various values of 7. 2 Also plotted 
in Figure [1] are points from a finite difference numeri- 
cal integration of equations (jTTj) - (114)1 . The code used 
to obtain these resu lts is a one-dimensional version of 
the ZEUS algorithm ([Stone fc NormarJ1992tlStone et all 

2 Comparison to the exact solution indicates that expression l|39| l 
is accurate to within 0.8% for 7 > 1.1. 
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119921 iTurner fe Stond l200ll) without diffusion. The ini- 
tial temperature was set to give (3q = cq/c = 10~ 4 
(this along with the value for ao sets the initial den- 
sity), and the piston velocity was set to co/(T + 1). 
Power-law f i ts to the Rosseland mean opacity were used 
(jBell fc Linlll994h . and the computational domain was 
set to L — , where ko(po,Tq) is the initial opacity. 3 
The slope was measured from a least squares fit to the 
velocity profile as a function of x after the wave front 
had propagated across the computational domain. 

Analytical profiles of velocity, density, temperature 
and a as a function of £ for 7 = 5/3 are shown in Fig- 
ures EKSl The fact that a ~ const, for ao 1 is consis- 
tent with expression (|19|) . which gives p oc T 3 for a» 1, 
i.e., a oc T 3 /p ~ const; it is simply the conservation of 
entropy in the high energy density limit. When 7 = 4/3, 
Ti = 4/3, c a = c m yjl + ao and expressions (|39|) and (|40| 
are exact. 

4.1. Isothermal Limit 

The considerations that lead to the approximate ex- 
pressions (|3"9")l and ([40]) break down when T\ ~ 1. In 
that case, and for a precise code comparison, the exact 
solution must be used (see the Appendix for details on 
its calculation) . The fluid quantity that is most sensitive 
to the breakdown of the approximate solution is c a . Fig- 
ures IH] and [7J show profiles of the adiabatic sound speed 
for 7 = 1.1 and 7 = 1.01. As the adiabatic index of the 
material approaches unity, the speed of sound in the fluid 
increases as it rarefies. This is generally associated with 
a value for Ti less than unity, 4 which from expression 
(|3Tj) occurs for a less than 



0.25 




-7-I' 
(44) 

where the latter expression is valid for 7-I < 1, This is 
equivalent to a ratio of energy densities near unity. Real 
solutions to the above expression can only be obtained 
for 7 < 13/12. The increase of c a in regions of rarefaction 
thus occurs for a material component with 



7-1 



> 24, 



(45) 



where / is the microscopic degrees of freedom. A more 
careful analysis (see the Appendix) confirms the / > 24 
result, with a slight modification to the value for a cr it 
(expression |A9| ) . The calculation in the Appendix also 
demonstrates that 1) the peak in c a seen in Figure [7J 
continues to grow with decreasing 7, with c a ,max ~ 
(7 — l) -1 / 2 (expression |A12j ). 2) the sound speed al- 
ways decreases for a > a cr u and 3) the portion of the 
solution that approaches the isothermal limit only exists 
for a < 0.714. 

3 This corresponds to an optical depth of unity across the com- 
putational domain. Without diffusion, the only constraint on the 
overall length scale is to = Lna 3> Po so that the coupling time 
scale is much less than the advection time scale. 

4 This demonstrates another breakdown in the analogy between 
a radiating fluid and a material fluid with 7 = Ti. If the latter 
were strictly true, Ti < 1 would imply a negative pressure; as it 
is, despite the somewhat unusual behavior of the sound speed, the 
pressure remains perfectly well behaved. 



Part of the reason for the unusual behavior of the sound 
speed is the increase of a as the fluid is rarefied (recall 
that this occurs for any 7 < 4/3). Figure [5] shows the 
profile of a for 7 = 1.01; even for ao < 1, the pres- 
sure ratio increases rapidly in the rarefaction region and 
the fluid quickly becomes dominated by radiation pres- 
sure, independent of the initial pressure ratio. Figure [§] 
shows the ratio of energy densities for the same set of 
parameters; this remains close to unity for a wider range 
of parameters. While the fluid near the quiescent state 
asymptotes to the hydrodynamic solution {T\ ~ 1), the 
fluid eventually (for a sufficiently large piston velocity) 
approaches the radiation-dominated solution (ri = 4/3) 
in which both radiation pressure and radiation energy 
density dominate. The increase of c a occurs in the tran- 
sition between these two asymptotic regimes where the 
fluid is dominated by radiation pressure but the mate- 
rial and radiation energy densities are comparable. One 
clear implication of these results is that RHD calcula- 
tions with 7 ~ 1 can produce drastically different results 
from a strictly isothermal calculation. 

4.2. Effects of Diffusion 

Figure [10] shows the velocity profile from finite differ- 
ence numerical integrations both with and without diffu- 
sion, for various values of /3qTq. The parameters for these 
calculations are 7 = 1.6667, j3o = 10 -4 and ao = 1. For 
sufficiently large values of /3o T 0: the diffusion is only mod- 
ifying the transition between the self-similar and con- 
stant portions of the fluid, as expected. 

5. CHARACTERISTICS AND REIMANN INVARIANTS 

As discussed in §101 of lLandau fc Lifshitd (|1987l ). all 
that is required to extend the results of the previous sec- 
tion to more general boundary conditions is the ability 
to express all of the fluid quantities as functions of one 
another (P = P[p], v = v[P], etc.). This is ensured for 
isentropic boundary conditions, and the generalization of 
expression (|35|) is 

x = t(v±c a ) + f(v), (46) 

where f(v) is an arbitrary function of the velocity. These 
solutions are referred to as simple waves, since the flow 
variables arc all functions of x ± c a t. The centered rar- 
efaction wave is a simple wave with f(v) = 0. 
The Riemann invariants for the combined fluid are 

J ± =v±f—. (47) 

J PCa 

They are conserved along the characteristic curves C+ = 
v±c a . As discussed in §104 of lLandau fc Lifshitd (|1987h . 
they are strictly conserved only for isentropic flow. For 
adiabatic flow, the perturbations Sv ± SP/ (pc a ) are con- 
served along the C± characteristics, and perturbations 
in the total entropy, Ss m + 5s r cx (SU + 5E — [7 17 + 
4E/3]6p/p)/U oc (l/[ 7 - 1] + 12a)<5T/T - (1 + 4a)Sp/p, 
are conserved along the characteristic Co = v. 

6. SUMMARY AND DISCUSSION 

This paper has regarded a radiating fluid in the dy- 
namic diffusion limit as ideal; diffusion is neglected for 
the bulk of the fluid since it only modifies regions of the 
flow in which gradients are large, in a manner analogous 
to other forms of dissipation such as viscosity and heat 
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conduction. A complete analysis has been performed 
for a centered rarefaction wave O, and expressions for 
the characteristic curves and Reimann invariants of ideal 
RHD have been provided (©. It has been shown that 
the much of the hydrodynamic analysis carries over by 
simply replacing the standard adiabatic quantities with 
the form that they take in a radiating fluid (@. The 
exact solution requires the numerical integration of ex- 
pression (|36[) and the implicit solution of equation (|35|) . 
The qualitative nature of the solution changes as the in- 
ternal degrees of freedom of the material become signifi- 
cant (/ > 24). In that case, if the initial energy content 
of the radiation is less the 90% of the internal energy of 
the material, the sound speed of the fluid increases as it 
rarefies. 

One application of 21 is as a test of numerical RHD 
codes; the solutions are nonlinear and the only term that 
has been neglected is the diffusion operator. One can 
ensure that the system of equations is in the dynamic 
diffusion limit simply by setting the length of the com- 
putational domain such that its optical depth satisfies 
r 3> c/c a . Other tests could be constructed by further 
extendi ng hydrodynamic results; to cite a couple of ex- 
amples, lLandau fc Lifshitzl <|1987f) include analyses of a 
uniformly accelerated piston and a centered rarefaction 
wave reflecting off a solid wall. The solution for a com- 
pressive wave before it steepens into a shock can be ob- 
tained from the negative branch of equations (|35[) and 

(25). 

Any application of the results of 21 to physical sys- 
tems must keep in mind the idealized nature of the anal- 



ysis. Since the density and temperature of the solutions 
change significantly from their quiescent values, there are 
likely to be large regions of parameter space for which the 
assumptions of dynamic diffusion break down. Mapping 
out the limits of their validity for realistic opacities would 
be a useful follow-on exercise to this work. For flows with 
/3t ~ 1, the modification of these solutions by diffusion 
in the transition regions could likely be obtained by an 
asymptotic analysis. 

The considerations of 21 could be used in the devel- 
opment of Godunov numerical schemes for RHD. An 
analysis of the hyperbolic nature of the full set o f RHD 
equations has been conducted bv iBalsara (1999a. b]); the 
radiation flux in these papers was regarded as a source 
term, however, so that only the material sound speed 
enters the expression for the C± characteristics. The 
analysis performed here points to a more self-consistent 
approach for numerically solving the equations of RHD 
in the dynamic diffusion limit. For calculations in which 
both advection and diffusion are important, one could 
still employ a Reimann solver for the hyperbolic portion 
of the equations while t reating the diffusion separately 
(iDragojlovic et alJl2006t ). 
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APPENDIX 
CALCULATION DETAILS 

The velocity can be calculated as a function of any of the other fluid variables by expressing the integral (|36)) in the 
form appropriate for that variable. The closed form expression of the integral in terms of a is given by 

1 + 12(7-1)™ , A1 , 

where c a is given by expression (|32[) . Numerical evaluation of this integral is straightforward for 7 > 4/3, but it 
becomes somewhat problematic when 7 < 4/3 due to the rapid increase of a. 5 Expressing the integral in terms of the 
logarithm of the density appears to be more robust, although other choices may be superior. Equation (f2"2"| can be 
inverted to express a in terms of the density: 



A — = W ( Aexp 
a 



A + (3 7 - 4) In £- 
Po. 



(A2) 



where W is the Lambert- W function (or product log) and 

A =12(7-1)00. (A3) 

The velocity is then given by 

v = / c a (rj) drj, (A4) 



where r\ = \n(p/po) and c a (rj) is obtained from expressions (|32| and (|A2|) . The flow profiles in the text were plotted 
with Mathematica, which provides a module for evaluating the product log, although one can generate its values 
based upon a simple recursion formula. With v(p) obtained numerically, equation ([33)) was solved with a root-finding 
algorithm to obtain p(£). All the fluid quantities as a function of £ then follow from the expressions given in the text. 

The unusual behavior of the sound speed as 7 1 discussed in the text can be obtained analytically as follows. 
Clearly an increase in the sound speed requires that there be one or more locations in the fluid for which dc a /d£ = 0. 
Differentiating expression (f3"2"]l . one finds that 

1 dc a (7 -!)(! + 4a) / 1+^4 5 - 3 7 + A \ 

c a da l + A V 4 (37-4)o 1 + A + (7 - 1)(1 + 4a) 2 / ' 1 ' 

5 To get the rarefaction solution, one must integrate between a = »o and a 5j ao when 7 Js= 4/3. 
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Setting this expression to zero gives an expression that is quadratic in 7 — 1 and cubic in a®. 

(1 + 12ea)(l + 12ea + e(l + 4a) 2 ) + 4(3e - l)a(2 - 3e + 12ea) = 0, (A6) 

where 

e = 7-l. (A7) 
The physical branch of expression (|A6|) regarded as an equation for e(a) is 

-1 - 68a + 32a 2 + (1 + 4a) V 1 + 224a + 448a 2 



48a(-l + 16a + 8a 2 ) 



(A8) 



The maximum of the above expression occurs at a = 0.714 and is 0.0862; i.e., 7 < 1.0862. This corresponds to 
/ > 23.2, or, restricti ng / to integer values, / > 24. 

Solving expression (|A6j) for a(e) yields two positive real roots when e < 0.0862. One of these corresponds to the 
peak value of c a and is approximately 

a cnt ~ ^^-2.22 + 1.15e + 0(e 2 ). (A9) 
e 



This corresponds to a ratio of energy densities 



— = Sea ~ 0.9. (A10) 



An estimate of the peak value of c a can be obtained by looking at its asymptotic form for e«1 and ae ~ 1: 



4<xJ- T^e- 2ae . (All) 

y 1 + 12ea v ' 



This has a maximum at ae = (1 + y/T)/12 = 0.304 given by 

0.307 



(A12) 



The other root varies between 0.125 < a < 0.714 for < e < 0.0862. This corresponds to a minimum in c a and 
marks the transition between the portion of the solution that approaches the isothermal limit and the portion in which 
the sound speed increases. For ao > a cr u, neither of the roots is accessible (since a > ato). For 0.714 < ao < ttcriti 
the first root is accessible but not the second. For ao < 0.714 both roots are accessible. 
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Fig. 1. — Velocity slope for 7 = 1.6667, 1.5, 1.4, 1.3333, 1.2, 1.1 and 1.01 (from bottom to top). The solid lines are the approximate 
analytical results and the points are results from a finite difference numerical integration with 512 grid points. 




Fig. 2. — Velocity profile for 7 = 5/3 and 00 = 0.03, 0.1, 0.3, 1, 3 and 10 (from top to bottom). Results outside of this range are close 
to the bracketing results. 
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Fig. 3. — Density profile for 7 = 5/3 and ciq = 0.01 (solid line), 0.3 (dotted line) and 10 (dashed line). Results outside of this range are 
close to the bracketing results. 




Fig. 4. — Temperature profile for 7 = 5/3 and ao = 0.01, 0.03, 0.1, 0.3, 1, 3, 10 and 30 (from bottom to top). Results outside of this 
range are close to the bracketing results. 




Fig. 5.— Profile of a for 7 = 5/3 and a = 0.01, 0.03, 0.1, 0.3, 1, 3, 10, 30, 100 and 1000 (from bottom to top). Results outside of this 
range are close to the bracketing results. 




Fig. 6. — Profile of c a for 7 = 1.1 and 00 = 0.001, 0.01, 0.1, 1, 10 and 100 (solid lines from left to right). Also shown are the 
hydrodynamic results for 7 = 1.01 (dotted line) and 7 = 4/3 (dashed line). 




Fig. 7.— Profile of c a for 7 = 1.01 and a = 0.001, 0.01, 0.1, 1, 3, 10, 30, 100 and 300 (solid lines from left to right). Also shown arc 
the hydrodynamic results for 7 = 1.01 (dotted line) and 7 = 4/3 (dashed line). 
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Fig. 8.— Profile of P r /P g for 7 = 1.01 and a = 0.001, 0.01, 0.1, 1, 3, 10, 30, 100 and 300 (from bottom to top). 




Fig. 9.— Profile of E/U for 7 = 1.01 and a = 0.001, 0.01, 0.1, 1, 3, 10, 30, 100 and 300 (from bottom to top). 
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Fig. 10. — Velocity profile from a finite difference numerical integration with 7 = 1.6667, ao = 1 and /3oro = 00 (solid line), 10 4 , 10 3 , 
10 2 and 10 (dashed lines in order of decreasing accuracy). For /3qto = 10, the temperature drop is sufficient to make fir ~ 1 at x ~ 0.1L. 



